CERN-TH/2003-42 



Resummed jet rates for e + e annihilation into massive quark^] 

Frank Kraus^j] and German Rodrigc0 

Theory Division, CERN 
CH-1211 Geneva 23, Switzerland 
(Dated: March 5, 2003) 

Expressions for Sudakov form factors for heavy quarks are presented. They are used to construct 
resummed jet rates for up to four jets in e + e~ annihilation. The coefficients of leading and next- 
to-leading logarithmic corrections, mandatory for a combination with higher order matrix elements, 
are evaluated up to second order in a s . 
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The formation of jets is the most prominent feature of perturbative QCD in e + e~ annihilation into hadrons. Jets 
can be visualized as large portions of hadronic energy or, equivalently, as a set of hadrons confined to an angular 
region in the detector. In the past, this qualitative definition was replaced by quantitatively precise schemes to define 
and measure jets, such as the cone algorithms of the Weinberg-Sterman [J type or clustering algorithms, e.g. the 
Jade 0j or the Durham scheme (k± scheme) : 4j. A refinement of the latter one is provided by the Cambridge 
algorithm [jj. Within the context of this paper, however, the Durham and the Cambridge algorithms are equivalent 
and they will be referred to as k± algorithm. A clustering according to the relative transverse momenta has a number 
of properties that minimize the effect of hadronization corrections and allow an exponentiation of leading (LL) and 
next-to-leading logarithms (NLL) stemming from soft and collinear emission of secondary partons. 

Equipped with a precise jet definition the determination of jet production cross sections and their intrinsic properties 
is one of the traditional tools to investigate the structure of the strong interaction and to deduce its fundamental 
parameters. In the past decade, precision measurements, especially in e + e~ annihilation, have established both the 
gauge group structure underlying QCD 0, H, El 0, El and the running of its coupling constant a s over a wide range 
of scales |l2j . In a similar way, also the quark masses should vary with the scale. 

A typical strategy to determine the mass of, say, the bottom-quark at the centre-of-mass (cm.) e nergy o f the 
collider is to compare the ratio of three-jet production cross sections for heavy and light quarks 0, Q, 15, 16, 171. 
At jet resolution scales below the mass of the quark, i.e. for gluons emitted by the quark with a relative transverse 
momentum k± smaller than the mass, the collinear divergences are regularized by the quark mass. In this region mass 
effects are enhanced by logarithms hi(mb/k±), increasing the significance of the measurement. Indeed, this leads to a 
multiscale problem since in this kinematical region also large logarithms ]n{y/s/k±) appear such that both logarithms 
need to be resummed simultaneously. A solution to a somewhat similar two-scale problem, namely for the average 
subjet multiplicities in two- and three-jet events in e + e~ annihilation was given in pa ]. 

These large logarithms can be deduced by inspection of the corresponding splitting functions for massive quarks [Tsl 
HH HE E3 an d of the boundaries for their integration over the energy fractions of the outgoing particles. The 
resulting integrated splitting functions exhibit the LL and NLL behaviour, and resummation is achieved by means 
of Sudakov form factors. Jet rates can then be expressed, up to NLL accuracy, via the integrated splitting functions 
and the corresponding Sudakov form factors obtained from them. Following the work of Catani et al. 4], the 
resummation of such logarithms will be discussed in this paper for the case of heavy quark production. Furthermore, 
the corresponding LL and NLL coefficients will be calculated. These coefficients are mandatory for the combination 
with next-to-leading order calculations of the three-jet rate, involving heavy quarks in e + e~ annihilation pH l23t Wk\ 
I25L l26l I27I l28l I29I ]. In fact, they exhibit the correct logarithmic behaviour and thus provide a good estimate for 
the size of mass effects. Such a matching procedure was first defined for event shapes in [3(j, and later applications 
include the matching of fixed-order and resummed expressions for the four-jet rate in e + e~ annihilation into massless 
quarks |3ll |32( • A similar scheme for the matching of tree-level matrix elements with resummed expressions in the 
framework of Monte Carlo event generators for e + e~ processes was suggested in [33| and extended to general collision 
types in j3^|. Finally, mass effects are briefly highlighted for the case of two- and three-jet events both for bottom 
quarks at LEP1 energies and for top quarks at a potential linear collider operating at cm. energies of 500 GeV. 
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To obtain expressions for the splitting function of the process a — > fe(pi) + c(p 2 ), which involves massive particles, 
we use the following Sudakov parametrization 
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where p M and n M are light-like vectors: p 2 = n 2 = 0, and q x is the space-like transverse momentum, p ■ q± = n ■ q± = 0, 
with q 2 = — q\ > 0. Furthermore, q — |q|. The quasi-collinear limit, as discussed in [l^, is obtained by rescaling q 
and to through q —>■ \q and to — ► Am, respectively. Taking the limit A — ► 0, and keeping only terms of order 1/A 2 in 
the squared matrix element we reproduce their result. It should be stressed, however, that this result is independent 
of the ratio q/m. In the limit discussed above, the squared amplitude at the tree-level fulfils a factorization formula, 
and it also contains naturally the two limits q <C to and q 3> to as particular cases. 
For the splitting process 



Q -> Q(pi) + ff(P2) , 

Q being a heavy quark, p\ = to 2 , the squared matrix element factorizes as 
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where the unregularized spin-averaged splitting function in D = 4 — 2e dimensions is given by 



(3) 



} (z,q) = C F 
= C F 



1 + z 2 
1-z 
1 + z 2 



-e(l-z)- 
-e(l-z)- 



Pl-Vi. 

2z(l — z)m 2 



1 — z ' q 2 + (1 — z) 2 m 2 

Analogously, for the g — > Q{p\) + Q(p 2 ) branching a similar factorization formula holds, where 
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As expected, these splitting functions match the massless splitting functions in the limit to — > for q 2 fixed. Finally, 
for the g —> g{p\) + g(p 2 ) branching 
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In the above equations, Cp — 4/3 and Ca = 3 are the structure constants of SU(3) in the fundamental and adjoint 
representation, respectively, and Tji = 1/2 is the normalization of its generators. For the purposes of this investigation, 
however, the splitting functions can be taken in D = 4 dimensions. 

Accounting for the number nj' h ^ of active light or heavy fermions (quarks), respectively, P 9 q can be replaced by 



Pgf( z i<l) = n f ■ PgQ( z :<l) 
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where the sum runs over all ri f flavours of heavy quarks. For massless particles the corresponding integrated splitting 



functions or branching probabilities T yield 
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The Sudakov form factors, which yield the probability for a parton experiencing no emission of a secondary parton 
between transverse momentum scales Q down to Q$, read 



Aq(Q,Qo) = exp 



A g {Q.,Qo) = exp 
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The fact that the above equations © and © are for massless particles only is reflected by the "propagator-like" 
structure 1/q. In contrast, for massive particles the propagator term is given by l/[q 2 + (1 — z) 2 m 2 ] for Q — ► Qg, 
and by l/[g 2 + m 2 ] for 5 — > instead. In order to compensate for this we define the integrated splitting functions 
involving heavy quarks through 
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Written in such a fashion, the Sudakov form factor for massive quarks is obtained by a mere replacement of the 
integrated splitting function according to Eq. (jSJ). Similarly, for every heavy quark Q occurring in a splitting g — > QQ, 
the integrated splitting function is modified correspondingly. In accordance with these replacements the running 
coupling constant a s changes as well, since the number of active quarks changes when passing a heavy quark threshold. 

The integrated splitting functions and the corresponding Sudakov form factors can be employed to calculate multi- 
jet rates in e + e~ annihilation in the kj_ schemes. In both schemes, the jet resolution parameter y^j is given by 
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where s = Q 2 is the cm. energy squared of the colliding electrons. For massive particles, it might be more suitable 
to replace the energies with the absolute values of the three-momenta in order to identify the relative transverse 
momentum properly, i.e. 
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However for our discussion this difference is subleading only. 

Rates for two-, three-, and four-jet events can be expressed by the integrated splitting functions and the Sudakov 
form factors as 
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where Qo now plays the role of the jet resolution scale in the k± algorithm, Qq = y cu tQ 2 , and Q is the cm. energy of 
the colliding e + e~. Single-flavour jet rates in Eq. I|13|l are defined from the flavour of the primary vertex, i.e. events 
with gluon splitting into heavy quarks where the gluon has been emitted off primary light quarks are not included in 
the heavy jet rates but would be considered in the jet rates for light quarks. 

In order to catch which kind of logarithmic corrections are resummed with these expressions it is illustrative to 
study the above formulae in the kinematical regime such that Q ^ m ^ Qq. Expanding in powers of a s , jet rates 
can be formally expressed as 
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where the coefficients cjy are polynomials of order I in L y = ln(l/y cut ) and L m = ln(m 2 / 'Qq) . At the given NLL 
accuracy it is sufficient to treat the running of a s to first order (one-loop), 



a s (q) = 
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l + /3 n ^ln(q 2 /Q 2 ) 
where the /^-function (3 n for n active quarks is given by 
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Alternatively, instead of using the transverse momentum q as the scale for a s in the Sudakov form factors, the 
expression showing up in the propagator terms, i.e. q 2 + (1 — z) 2 m 2 for the heavy quark splitting Q — > Qg, or q 2 + m 2 
for the gluon splitting into two heavy quarks g — > QQ, can be chosen. However, this is a subleading effect. 
The coefficients for the first order in a s are given by 
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all coefficients for higher jet multiplicities being 0. For second order a s with n active flavours at the high scale, the 
LL and NLL coefficients read 



.(2) _ 



„(3) _ 
-24 ~ 



-.CI I / Li 



"24 



= ^Ci L 2 - L; 



— C f Ca ( - L 4 m 



lo 9 10 (yJ 



lo 9w(y C ut) 



FIG. 1: Effect of a 6-mass of 5 GeV on the single-flavour two- and three-jet rate at LEP1 energies as a function of the jet 
resolution parameter in the Durham scheme. In the left plot this effect is treated through the full inclusion of masses into the 
splitting function, see Eq. I1LH . whereas in the plot on the right hand side this effect is modeled through the dead cone, see 
Eq. CSS). 
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Terms ~ (f} n — f3 n -i) in the NLL coefficients are due to the combined effect of the gluon splitting into massive quarks 
and of the running of a s below the threshold of the heavy quarks, with a corresponding change in the number of 
active flavours. With our definition of jet rates with primary quarks the jet rates add up to one at NLL accuracy. 
This statement is obviously realized in the result above order by order in a s . 

The corresponding massless result |4( is obtained from Eqs. I|17|l and l|18fl by setting L m — > 0. Notice that Eqs. 117|) 
and l|18l) are valid only for m ^> Qa and therefore m — > does not reproduce the correct limit, which has to be smooth 
as given by Ea. H13|) . Let us also mention that for Q > m there is a strong cancellation of leading logarithms and 
therefore subleading effects become more pronounced. 

An approximate way of including mass effects in massless calculations, that is sometimes used, is the "dead cone" 
approximation. The dead cone relies on the observation that, at leading logarithmic order, there is no radiation of 
soft and collinear gluons off heavy quarks. This effect can be easily understood from the splitting function in Eq. fijl. 
For q <C (1 — z) m the splitting function is not any more enhanced at z — ► 1. This can be expressed via the modified 
integrated splitting function 

Iq c (Q, q, m) = Tq(Q, q, m = 0) + 2C F In (£) 9(m - q) . (19) 

To obtain this result the massless splitting function has been used, which is integrated with the additional constraint 
z > 1 — q/m. We also compare our results with this approximation. 

The impact of mass effects can be highlighted by two examples, namely by the effect of the 6-mass in e + e~ 
annihilation at the Z-pole and by the effect of the i-mass at a potential linear collider operating at a cm. energy of 
500 GeV. With m b = 5 GeV, M z = 91.2 GeV, and a s (M z ) = 0.118, the effect of the 6-mass at the Z-pole on the two- 
and three-jet rates is depicted in Fig. ^ Clearly, using the full massive splitting function, the onset of mass effects in 
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FIG. 2: The effect of a t-mass of 175 GeV on the single-flavour two- and three-jet rate as a function of the jet resolution 
parameter, at a potential linear collider operating at cm. energies of 500 GeV. Again, in the left plot this effect is treated 
through the full inclusion of masses into the splitting function, see Eq. (I1U1 . whereas in the plot on the right hand side this 
effect is modeled through the dead cone, see Eq. 11911 . 

the jet rates is not abrupt as in the dead cone case and becomes visible much earlier. Already at the rather modest 
value of the jet resolution parameters of y cu % — 0.004, the two-jet rate, including mass effects, is enhanced by roughly 
4% with respect to the massless case, whereas the three-jet rate is decreased by roughly 3.5%. For even smaller jet 
resolution parameters, the two-jet rate experiences an increasing enhancement, whereas the massive three-jet rate 
starts being larger than the massless one at values of the jet resolution parameters of the order of y cu t ~ 0.001. The 
curves have been obtained by numerical integration of Eq. (|13|l . Furthermore, in order to obtain physical result the 
branching probabilities have been set to one whenever they exceed one or to zero whenever they become negative. 

While in the case of bottom quarks at LEP1 energies the overall effect of the quark mass is at the few-per-cent 
level, this effect becomes tremendous for top quarks at a potential linear collider operating at 500 GeV; see Fig. [21 
Owing to the size of mass effects, the difference between the LL treatment through the dead cone and the full NLL 
result becomes visible, reflecting the fact that Q ~ m and the respective logarithms canceling each other. In other 
words, for this case a full fixed-order treatment is mandatory. 

In this paper Sudakov form factors involving heavy quarks have been employed to estimate the size of their mass 
effects in jet rates in e + e~ annihilation into hadrons. These effects are sizeable and therefore observable in the 
experimentally relevant region. In addition, the coefficients for the leading logarithms, both in the jet resolution 
parameter and in the quark mass, have been deduced up to second order in a s . They are mandatory for the 
combination with fixed-order calculations of the two-, three- and four-jet rates, and resummed expressions as obtained 
by Sudakov form factors. The matching between fixed-order calculations and resummed results will be presented in 
a forthcoming article. 
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